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ABSTRACT 

Wc describe the motivation, field locations and stellar selection for the ARGOS spec- 
troscopic survey of 28,000 stars in the bulge and inner disk of the Milky Way galaxy 
across latitudes of b = —5° to —10°. The primary goal of this survey is to constrain the 
formation processes of the bulge and establish whether it is predominantly a merger 
or instability remnant. From the spectra (R = 11,000), we have measured radial ve- 
locities and determined stellar parameters, including metallicities and [a/Fe] ratios. 
Distances were estimated from the derived stellar parameters and about 14,000 stars 
arc red giants within 3.5 kpc of the Galactic centre. In this paper, we present the 
observations and analysis methods. Subsequent papers (III and IV) will discuss the 
stellar metallicity distribution and kinematics of the Galactic bulge and inner disk, 
and the implications for the formation of the bulge. 

Key words: Galactic bulge, rotation profile, [Fe/H]. 



1 INTRODUCTION 

The Abundances and Radial velocity Galactic Origins Sur- 
vey (ARGOS) is a large spectroscopic survey of about 28,000 
stars in the Galactic bulge and surrounding disk, using the 
AAOmega fibre spectrometer on the Anglo- Australian Tele- 
scope. With this survey, we aim to investigate the metallic- 
ities, kinematics and density distributions of the stars in 
the inner region and use this information to understand 
the bulge. Near-infrared images of the Galactic bulge (e.g. 
the COBE/DIRBE images) show its peanut /boxy structure. 
Peanut /boxy bulges are widely believed to be generated via 
bar-forming and bar-buckling instabilities of the disk, and 
the main motivation of this survey is to understand how 
the Galactic bulge formed. Did it form primarily through 
disk instability or did mergers play a significant role in its 

formation? 

Several authors (e.g. lAtha nassoula 2005: 



Martinez- Vafpuesta. Shlosman. fe Heller 20061 ; Shen et all 
201Ch have made N-body simulations of peanut/boxy bulge 



formation through bar-forming and bar-buckling instability 
of the disk. One of the goals of the ARGOS survey is 
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to see whether the existing N-body models provide a 
good representation of the bulge properties, in terms of 
the observed kinematics and peanut /boxy structure. The 
ARGOS data provide 3D positions and radial velocities for 
the stars, and also their [Fe/H] and [a/Fe] abundances. The 
combination of positional, kinematical and chemical data 
will be useful for testing future models of bulge formation 
that include chemical evolution. 

A generic prediction of the N-body models is that 
instability-generated bulges are expected to show cylindri- 
cal rotation, unlike the rotation field predicted for the more 
spheroidal classical bulges. The prediction of cylindrical ro- 
tation for the boxy bar/bulge is easily checked from radial 
velocity surveys of the bulge, and has already b een con- 
firmed by the BRAVA group (|Howard et al.l 120091 ). Testing 
for the prese nce of a classical component of the bulge is 
more difficult. ISaha. Martinez- Vafpuesta. fc Gerhard! l|2012h 
showed that a small classical bulge in the presence of a boxy 
bar/bulge would itself be spun up into cylindrical rotation 
and would therefore lose its distinct kinematical signature. 

Recently, a split red clump in the Galactic bulge has 
been discovered (see ISaito et all 1201 ll ; iNataf et alj 20101 ; 
iMc William fc Zoccalill2010h . The split is seen as a bimodal 
magnitude distribution of the red clump stars at higher lati- 
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tudes near the minor axis of the bulge, and is associated with 
the boxy/peanut structure of the bulge. In projection, the 
boxy/peanut bulge shows an underlying X-structure. Our 
data are well suited to investigating the stellar populations 
present in the X-structure. The first pap er of this series i s 
on the X-structure observed in the bulge (|Ness et alj|2012ft . 
and showed that this characteristic structure of the bulge is 
defined by stars with metallicities [Fe/H] > —0.5. The more 
metal-poor stars do not appear to show the split red clump 
and are therefore probably not involved in the boxy/peanut 
structure. 

If the Galactic bulge formed from instability of the inner 
disk many Gyr ago, then the bulge would be made up mainly 
of early thin and thick disk stars that have been trapped 
within the bulge structure. These stars would present as 
a fossil of the early inner Galaxy, and we could expect 
to see evidence of these trapped early components in the 
metallicity distribution and kinematics of the bulge. The 
relative contributions of the early components are likely 
to vary across the bulge, and this variation may give in- 
sight into the mapping of the early disk into the bulge 
via the instability. We therefore chose to survey stars in 
the bulge and also in the surrounding thin and thick disk, 
to see whether the metallicity distribution across the in- 
ner regions of the Galaxy is consistent with the instabil- 
ity picture. In this context, the Galact ic bulge appears to 
show a vertical abundance gradient (e.g. Minniti et al.|[l995l ; 
IZoccali et all 120081 ; iBabusiaux et al.ll201Ct l" which could be 
interpreted as evidence for a dissipative formation. It is of in- 
terest to investigate this abundance gradient in more detail. 
Futhermore, the [a/F el properties of the bulge have been 
much discussed (e g iFulbright. McWilliam. fc Rich! 120071; 
Zoccali et alll2007l; iMelendez et al.ll2008l ; lAlves-Brito et all 



20ld ; I Johnson et alj|201lh . with indications from some au 



thors that ev en the most metal-r ich stars in the bulge are 
a— enhanced. iBensbv et ail ||201Ch also find evidence for en- 
hanced [a/Fe] in a limited sample of inner disk giants (typi- 
cally > 1 kpc above the plane). The [a/Fe] ratios for stars in 
the inner Galaxy are an important diagnostic of the rate of 
star formation in this region. We therefore wanted to mea- 
sure the [a/Fe] ratios for the large sample of ARGOS bulge 
stars. 

There is a possibility that the first stars are concen- 
trated to the inner region s of the Galaxy, as indicated in cos- 
mological simulati ons (e.g lDiemand. Madau. fc Moorell2005l ; 
iBrook et all 12007ft . These stars are of great interest for un- 
derstanding the early chemical evolution of the Galactic en- 
vironment. The selection criteria for the survey have been 
defined to ensure we did not exclude such metal-poor giants. 

Along lines of sight passing through the inner re- 
gion of the Galaxy are stars from all the major Galactic 
components: thin and thick disk, bulge and stellar halo. 
Some of these are likely to be involved in the peanut/boxy 
bulge structure and some are not. Contamination by fore- 
ground and background disk stars is an issue for surveys 
of stars in the Galac tic bulge. Galactic models like Galaxia 
l|Sharma et al.ll201lft show that samples of red giants in the 
direction of the bulge can include significant fractions of 
stars of the thin and thick disk lying in front of and behind 
the bulge. Stars of the inner disk are of interest for under- 
standing how the bulge relates to its environment, particu- 
larly if it did form via disk instability. For specific studies 



of metallicities and kinematics in the bulge, however, we 
want to choose stars such as red clump giants for which 
accurate (~ 10%) distances can be estimated. With these 
distances, we can then determine which stars lie within the 
bulge and identify the stars which are in the foreground and 
background and derive a relatively uncontaminated picture 
of the distribution of kinematics and chemical properties 
throughout the bulge. For example, with the selection cri- 
teria described below, designed to select mainly red clump 
giants along sight lines through the bulge, we found that 
about 2,500 of our 28,000 stars are foreground dwarfs and 
about 14,000 stars lie within the bulge region (Galactocen- 
tric cylindrical radius Re < 3.5 kpc). The remaining 11,500 
stars are in the foreground or background and outside the 
inner region. 

The overall outcome of the survey is that the bulge ap- 
pears to be consistent with formation via bar forming and 
bar buckling instability of the early disk. The dominant com- 
ponent of the peanut/boxy bulge is the early old thin disk, 
with a clear contribution from its more metal-rich and kine- 
matically colder component, similar to the younger, more 
metal-rich and colder component seen in the thi n disk near 
the Sun at the present time (e.g. iHavwoodl 20081). 

The first paper of this series (|Ness et al.ll2012ft showed 
that the peanut/boxy structure of the bulge is defined by 
stars with metallicities [Fe/H] > —0.5, and its structure and 
kinematics are well represented by N-body models of bulge 
formation via disk instability. This paper, ARGOS II, deals 
with the observational methods and strategy for the ARGOS 
survey. Paper III will present the components of the metal- 
licity distribution for the bulge, and paper IV presents the 
kinematics of its sub-populations. Paper V will discuss the 
mapping of the early disk into the bulge by the instability. 



2 THE DATA 

In this section, we describe the selection criteria for the stars, 
and provide details of the AAOmega observations. 



2.1 Sample Selection 

Our goal is to select red giants in a grid of fields over 
the bulge and inner disk, with an apparent magnitude and 
colour range that includes clump giant stars all along the 
line of sight through the bulge, from the near side to the 
far side, and surrounding inner disk. The distance range is 
about 4.5 to 13 kpc from the Sun (see Figure [T]). We also 
wish to include metal-poor giants in the inner Galaxy, as de- 
scribed above, because of their interest as possible first stars. 
With the above goals and criteria in mind, we minimise the 
numbers of the brightest and coolest giants for which abun- 
dances are difficult to measure reliably from optical spectra. 
Foreground bluer main sequence stars should be excluded 
without excluding metal-poor giants along the line of sight. 

In order to complete the survey in a realistic time, we 
used the sparse sampling shown in Figure [2] This grid of 
field centers covers the southern bulge between b = —5° 
to —10°, extending out in longitude into the surrounding 
thin disk and thick disk. For the stars in the bulge re- 
gion, we choose field centres at latitudes \b\ ^ 5°, because 
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Figure 1. This sketch in the xy-planc illustrates the range of 
distance from the Sun that we wish to cover in each field with 
our selection of clump giants. 
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Figure 2. Galactic latitude and longitude for the 28 2° fields in 
our survey. The rectangle shows the approximate extent of the 
boxy bulge. The filled circles are fields for which about 1000 stars 
were observed. The open circle is our one incompletely observed 
field (~ 600 stars). 



the mean reddening increases rapidly closer to the Galac- 
tic plane. The field centres also include two low-latitude 
fields in the disk at larger longitudes in windows of low 
reddening, and three relatively low extinction fields in the 
northern bulge for comparison with the southern fields as a 
check on north-south symmetry. About 1000 stars were ob- 
served in each field. Stars were selecte d from the Two Micron 
All Sky Survey (2MASS) catalogue (|Skrutskie et al.l 120061 1 
over the magnitude range K = 11.5 to 14.0, with J,K er- 
rors < 0.06, and flags set to reduce blends, contamination, 
stars of lower photometric quality and non-stellar objects 
(Qflg=AAA, Bflg=lll, Cflg=000, Xflg=0, Aflg=0). 

The red clump giants have a well-defined absolut e mag- 
nitud e distribution, with a mean M Ko = -T61 ± 0.03 (|Alvesl 
l2000l l and a standard deviation of 0.22 mag. This allows a 
relatively accurate estimate of their distances (standard er- 
ror about 0.9 kpc at a distance of 8 kpc). The clump gi- 
ants lie in the plume near (J — K) — 0.65 in Figure 
With the range of apparent magnitudes selected for our sur- 
vey, the clump stars are distributed over the whole line 
of sight through the bulge (Figure [TJ. Some intrinsically 
brighter/fainter giants on the far/near side of the bulge 
are included in the sample, although the magnitude cut at 
K = 11.5 excludes most of the brighter and cooler giants 
for which it is more difficult to measure abundances. A blue 
colour cut was made to exclude the plume of bluer fore- 
ground stars seen in Figure The colour of this cut is at 
(J — K) = 0.38 so depends on the mean reddening in the 
field. The stars in the sample are then mainly clump gi- 
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Figure 3. Colour Magnitude Diagram from 2MASS for the bulge 
field at I = —10°, b = —5°. The red giant clump at the distance 
of the Galactic centre is near (J — K) a = 0.65, K = 13.2. The 
broken lines show the region from which our stars are selected. 
This region includes mainly red clump giants along the line of 
sight through the bulge, plus some brighter and fainter giants 
lying mostly on the far and near side of the bulge respectively. 
The plume of foreground disk stars to the left of the CMD is 
excluded by the colour cut in (J — K) but metal-poor giants 
(i.e. [Fc/H] < -2.0) are included as shown. 



ants, plus a contribution from brighter and fainter giants 
and some metal-poor giants which come mainly from the 
inner regions of the metal-poor halo of the Galaxy. We can 
expect some contamination from nearby G and K dwarfs, 
but they can be easily removed using their inferred stellar 
parameters. 

The stellar parameters are derived from our pri- 
mary wavelength region around the Ca-triplet. The high- 
est spectral resolu tion available with the AAOmega system 
(|Sharp et alj|2006t l at the Ca-triplet region (R = 11,000) is 
well suited to our goals: it provides accurate radial veloci- 
ties and abundances for a few individual elements including 
[Fe/H] and a-elements (see Figure [3J|. The stars are suffi- 
ciently bright to observe in b r ight t ime. 

The mean ISchlegel et ail (|l99St ) reddening for each AR- 
GOS field was adopted to make the blue cut. The reddening 
ranges from E(J - K) = 0.06 to 0.17 in the b = -10° fields, 
0.12 to 0.24 at b = -7.5°, 0.24 to 0.39 at b = -5° and 0.12 
to 0.21 in the northern fields and disk well away from the 
bul ge. Most of our star s lie at ^-heights above 300 pc, so 
the lSchlegel et all (|l99&f ) red dening is likely t o be a realistic 
estimate. We have used the ISchlegel et al.l (|l998l ) redden- 
ing throughout this s tudy, which was completed before the 
iGonzalez et alj (12011) reddening ma p for the bulge became 
available. The lGonzalez et alj (|201ll ) reddening in our field 
at (l,b) = (0° -5°) i s about 0.1 lower in E(J - K) than the 
ISchlegel et alj (|l998r ) reddening. The observed colour of the 
low density region in the CMD at (J—K) = 0.38, where the 
foreground star plume and the clump star sequence meet, 
provides a check on the mean reddening in the field. In 
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Figure 4. A typical AAOmcga stellar spectrum of a red giant 
star in the Ca-triplet region after continuum normalisation. Lines 
of some elements are indicated. All species are singularly ionised 
except Ca. The star shown here has stellar parameters of T e ff = 
4600 K, logg = 2.8, [Fe/H] = and a S/N = 70 per resolution 
clement. 

all fie lds, its colour was consistent with the ISchlegel et al.l 
(| 19981 ) reddening. The blue cut is sufficiently blue such that 
metal-poor giants in the field are not excluded. 

Towards the faint end of the selection interval in K- 
magnitude, the 2MASS catalog in the bulge fields becomes 
rapidly incomplete. We wanted to have clump stars all 
along the line of sight through the bulge, so chose roughly 
equal numbers of stars per unit magnitude in each field. 
It was also desirable to have a relatively small magnitude 
range for each fibre setup to reduce the effects of fibre- 
to-fibre scattering. Because the Ca-triplet region is close 
to the photometric /-band, we estimated the I-magnitude 
for each star from its JK photometry, using the transfor- 
mation I = Jo + 1.095 ( J - K ) + 0.421 E(B - V), de- 
rived from iBessell &l Brett! (| 19881 ) and adopting the ratios 
E(J - K)/E(B — V) = 0.535, A K = 0. 347 E(B - V) 
and Ai/ Ay — 0.601 and reddening values |Schlegel et al.l 
119981 ) taken from the NASA/IPAC Extragalactic Database 
(NED). Each field typically had many thousands of stars 
satisfying the selection criteria in Figure [3] 

2.2 Observations 

The spectra were acquired with the fibre-fed AAOmega sys- 
tem on the Anglo- Australian Telescope in May-August of 
2008 to 2011. AAOmega is a dual-beam spectrometer with 
about 400 fibres that can be positioned within a 2— degree 
diameter field. With AAOmega (allowing for broken fibres), 
we can observe about 350 stars at once, plus 25 sky posi- 
tions and guide fibres. The red beam provides our primary 
spectral region, using the 1700D grating with a resolution 
of about 11,000 in the Ca-triplet region 8400 - 8850A. This 
region is relatively free from telluric lines and contains spec- 
tral lines of Fe I, Ti I, Ca II, Al I, Mg I, Si I and O I for 
abundance measurements: see Figure 2] Of the optical re- 
gions of the spectrum that we could have chosen for this 
project, the Ca-triplet region is least affected by moonlight 
and it was possible to conduct the program during bright 
lunar phase. On the other hand, the region is not ideally 
suited to estimating the spectroscopic stellar surface gravi- 



ties which are needed to determine which stars are clump gi- 
ants and to estimate the distances for red giants that are not 
on the clump. The typical signal to noise ratio (S/N) for our 
Ca-triplet spectra is about 75 per resolution element. The 
blue beam was used to acquire low resolution spectra (R ~ 
3000) with the 1500V grating in the region 5000 - 5600A. 
The main purpose of these low resolution spectra is to re- 
ject foreground dwarf stars from our sample of more distant 
giants, using the Mgb/MgH feature. 

The spectra from the fibres in each of the red and blue 
beams are imaged on to a 2K x 2K E2V CCD. For each 
field, 25 fibres were used to measure the sky spectrum. Wave- 
length calibration spectra were acquired before and after the 
exposures for each fibre setup, and flat field spectra were 
taken for each setup in order to trace the spectra of the 
stars and to remove small fringing effects and pixel-to-pixel 
sensitivity variations. 

Figure [2] shows the locations of the 28 fields selected for 
this program. Spectra were acquired for about 28,000 stars, 
using 42 nights in the winters of 2008 — 2011. Our goal was 
to acquire spectra of about 1000 stars in each field. This was 
achieved, except for a single field at (I, b) — (20°, —7.5°) for 
which we have about 600 stars. 

The /—magnitudes of the selected stars were mostly 
between about 13 and 16 and were used to allocate each 
star to one of three one-magnitude intervals for observa- 
tion. These magnitude intervals were approximately from 
/ = 13 — 14, 14 — 15 and 15 — 16, depending on the mean 
reddening in the individual fields. In each field, we selected 
about 1000 stars, uniformly distributed over the three mag- 
nitude intervals and uniformly distributed over the 2° field 
so that every available fibre could be allocated to a star. 
The typical total exposure times for the three magnitude 
intervals were (60, 120, 150) minutes respectively, acquired 
in several sub-exposures, so the S/N is usually higher for the 
brighter stars. The clump stars in these /—band intervals are 
at mean distances of about (5.5, 8.5, 10.5) kpc respectively, 
covering the entire line of sight through the bar/bulge at 
longitudes between +20 and —20 degrees. For each field, 
25 sky positions, which are free of visible stars on the red 
Digital Sky Survey survey images, were selected visually. 

The S/N achieved in each field for our primary spectrum 
over the Ca-triplet region is given in Table [T] where the 
different magnitude ranges are indicated using 1, 2, 3, from 
the brightest to the faintest. The code for the field names is 
the longitude and then the latitude, where the m and p are 
for minus and plus values of I and 6: e.g. m0ml0_3 denotes 
the fibre setup at (l,b) = (0°,— 10°) and at the faintest 
magnitude range. 



3 DATA REDUCTION 

The first stage of the reduction used the standard 2dfdr 
pipelin^E This program subtracts bias and scattered light, 
extracts the individual stellar spectra, corrects for the flat 
field, calibrates wavelength and fibre throughput and sub- 
tracts sky. We selected the optimal sky subtraction option 
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Table 1. Mean S/N per resolution element in the Ca-triplet re- 
gion for each field. 
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which gives much better sky subtraction, and used a 4th 
order polynomial for fitting the arc lamp lines. 

The automated sky subtraction did not work well for 
low S/N data or on nights that were partially cloudy. There- 
fore, the sky subtraction was done externally with a python 
program that subtracted the sky on a star by star basis us- 
ing the sky emission lines for normalizing the throughputs 
of the star and sky fibres. The continuum sky was measured 
from the dedicated sky fibres; this improved the sky sub- 
traction and reduced the internal error of the radial velocity 
measurements via cross correlation with synthetic template 
spectra. For estimating the stellar parameters, the wave- 
length intervals of the strong sky lines were excluded from 
the analysis to ensure that the estimates were not affected 
by errors in sky subtraction. This turned out to be an un- 
necessary precaution, because the parameter estimates were 
unchanged whether or not the regions of the strong sky lines 
were included. The derived stellar parameters were similar 
for spectra reduced entirely with 2dfdr and spectra for which 
the external sky subtraction was used, except for the effect 
on the error of the radial velocities mentioned above. 

Scattered light removal is critical to ensure the correct 
continuum level. The scattered light removal with 2dfdr was 
checked using independent IRAF (dofiber) reductions. The 
IRAF results were consistent with those from 2dfdr. 

After data reduction, the spectra were normalized by 
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Figure 5. Spectrum of a representative star from Figure 2 show- 
ing the c ontinuum regions i n the blue shaded areas, identified 
using the lHinkle et aL [ [120051 ) atlas, used for normalisation of the 
AAOmega spectra 



fitting either a low order spline or Legendre polynomial to 
the median values of continuum regions identified from the 
Hinkle Arcturus Atlas (see Figure [5}. A small fraction of the 
spectra were not well fit using a single spline function; for 
these spectra we tried both functions (a 4th order Legendre 
polynomial and a 3rd order spline) and the better one was 
adopted. 



4 ANALYSIS 

All spectra were inspected by eye for an initial culling of the 
data to remove spectra that are significantly damaged by 
bad pixels or cosmic rays in critical regions and to remove 
double lined binaries. This reduced the original sample to 
about 26,000 stars which were passed through the x 2 pro- 
gram and analysed. A further 500 stars with S/N < 35 (per 
resolution element) were removed from the analysis, reduc- 
ing the final sample to 25,500 stars with a mean S/N of 
about 75 ± 20. The foreground cool dwarfs were eliminated 
from the sample using the x 2 program described below to 
estimate their gravities; this step follows a preliminary cut 
using the gravity-sensitive MgH and Mgb features in the 
blue region of the spectra to identify dwarfs. 



4.1 Overview: Measuring Stellar Parameters 

Stellar parameters (logg, [Fe/H] and [a/Fe]) were deter- 
mined using a minimum \ 2 program written for this dataset. 
The program compares each observed spectrum to a large 
grid of synthetic spectra, generated u sing the LTE stel- 
lar synthesis program MOOG 2010 dSnedenl Il973, and 
subsequent updates) from the ICastelli fc Cacciaril ( 200lh 
ID model atmospheres with no convective overshooting 
|Castelli fc Kuruca 120041 ') . The grid covers the parameter 
space T eff = 3750 - 6250 K, [Fe/H] = -5.0 to 1.0, logp 
= 0.0 to 5.0 and [a/Fe] = 0.0 to +1.0. The grid steps are the 
steps of the original opacity distribution functions (250 K in 
Teff, 0.5 in log ,g, 0.5 in [Fe/H] and 0.10 dex in [a/Fe]), for 
a microturbulent velocity of 2.0 km s _1 . An assumption of 
a constant microturbulent velocity may be a source of error 
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in our analysis. However, given we are observing mostly red 
clump stars which have a similar log g, any error introduced 
by this adopted value of 2.0 km s _1 should be systematic. 

The temperature is not estimated by the minimum x 2 
procedure because of the usual degeneracies involved in at- 
tempting to measure all of the stellar parameters from spec- 
tra in the Ca-triplet region. The temperature is determined 
from the (J — K) a colours of our giant sta rs using the cal- 
ibration by iBessell. Castelli. fc Plezl f|l998l 'l . The tempera- 
ture error, estimated to be ±100 K, comes from uncertain- 
ties in (i) the interstellar reddening, and (ii) the 2MASS 
apparent magnitude. All of our selected stars have JHK 
errors < 0.06 mag. Comparison with a set of calibration 
stars from the literature, presented in the Appendix, shows 
that our mean values of abundance and \ogg are close 
to those of the calibration stars. I f we had adopted the 
lAlonso. Arribas. fc Martinez- Rogerl (|l999h empirical tem- 
perature - colour calibration, the temperatures would be 
lower by about 100 K in the mean (up to 200 K for the most 
metal-poor stars and zero for the more metal-rich stars with 
[Fe/H] > 0.0 for which the temperatures from the two cali- 
brations agree). The mean change in temperature is within 
the errors of the calibration and would typically shift the \ 2 
[Fe/H] abundances down by about 0.1 dex in the mean and 
up to about 0.2 dex for the most metal-poor stars. 

The calibration stars include 16 bulge stars w i th [Fe /H] 
between 0.1 and —1.2 observed by IZoccali et al.l (|2008T ) at 
high resolution, and open and globular cluster stars ob- 
served by Da Cost a et al. (2011, private communication) and 
lLane et al.l (|201 lh with the same instrument and grating. 
The clusters include the open cluster Melotte 66 with [Fe/H] 
= -0.3, and 6 globular clusters: 47 Tuc ([Fe/H] = -0.75), 
NGC280 ([Fe/H] = -1.3), NGC1904 ([Fe/H] = -1.6), 
NGC6658 ([Fe/H] =-1.7), NGC2290 ([Fe/H] = -1.9) and 
NGC5034 ([Fe/H] = -2.1). The total range in [Fe/H] for 
the calibration stars is 0.1 to —2.10 . 

The solar and Arcturus spectra from iHinkle et al.l 

|2005l ) were also included in the calibration, after convo- 
lution to our resolution of 11,000. The errors in our parame- 
ters, as estimated from comparison of the \ 2 program results 
with the calibration star parameters, are ±100 K in T e //, 
±0.30 in logg, ±0.09 in [Fe/H] and ±0.10 in [a/Fe]. The 
comparison of the \ 2 an d literature results for these calibra- 
tion fields is given in the Appendix (see Figures |A 1 11 A21 |A3[) . 
For clarity, in Figure |A"T1 and Figure [A"2l the cluster data are 
represented by a single point with the observed dispersion 
but our errors in the measured parameters have been esti- 
mated from the entire calibration data set. The one-sigma 
errors for each parameter are indicated by the horizontal 
dashed lines. The metallicities of the globular clusters are 
adopted from the Harris database and the a-enhancement 
is a median compilation of literature sources for each clus- 
ter. For the internal uncertainties, a temperature change of 



^ References for cal ibration st ars: 47 Tuc: 

IMcWilliam fc Bernstein! teOOSf ). NGC 288: IShetrone fc Keand 
teOOCf). NGC 1904: iGratton fc Ortolanil Jl989f) ; 
Kirbv. Guhathakurta, fc Snedenl teOQgf), Melott e 66: 
Gratton fc Contarinil <1994h; ISestito et all <200Sh NGC 
2298: IMcWilliam. Geisler. fc Rich! dl992T) . Arcturus: 



± 100 K corresponds to a logg change of ± 0.3 dex and a 
[Fe/H] change of ± 0.1. 



iRarm'rez fc Allende Prietdl201irT 



4.2 Linelist 

The adopted linelis t comes from 

iKirbv. Guhathakurta. fc Snedenl J2008f). originally as- 
sembl ed from both t h e VA LD Kupka fc Rvabchikoval 
<|1999D and the iKuruc 

3 Ell) 

databases of atomic tran- 
sitions. The VALD database was used for a ll neutral and 
singly ionized atoms and the iKurucj |l992T i database for 
CN, C2, MgH as well as hyperfine transitions. The log gf 
values of the Fe lines and lines in surrounding regions were 
adjusted so that our synthetic spectra gave as good a match 
as possible to both the solar spectrum and the Arcturus 
spectrum (Hinkle Atlas) at our resolution of 11,000 and 
also to our spectra of the 16 bulge stars from Zoccali (see 
Appendix). The solar a bundances were take n from MOOG 
which uses all of the lAnders fc Grevessd jl989h values 
except for Fe. 

At this resolution there is some blending of lines, pre- 
dominantly with CN bands. The empirical log gf value ad- 
justments are specific to this resolution and reflect the level 
of blending. The synthetic spectra are not a perfect match 
to the observed spectra for any particular set of stellar pa- 
rameters; the log gf values were adjusted to give a good fit 
over the abundance range —1.2 < [Fe/H] < 0.1 covered by 
the Sun, Arcturus and the Zoccali stars. We then checked 
the stellar parameters returned by the \ 2 program for our 
more metal-poor globular cluster stars at lower [Fe/H]. (We 
recall that the [Fe/H] values come only from Fe lines.) No 
obvious deviations were seen (see Figure [A"2)l . We conclude 
that the derived stellar parameters for giants are on a con- 
sistent abundance scale in the range —2.1 < [Fe/H] < 0.1. 
Outside this abundance range, our stellar parameters are 
unconstrained by empirical calibration. 

The derived stellar parameters from our \ 2 program 
for the Sun are logg = 4.39 ± 0.30, [Fe/H] = -0.02 ± 0.09, 
[a/Fe] = ±0.04 ± 0.10 for a fixed T eff = 5780 K. For 
Arcturus the parameters are logg = 1.7 ± 0.30, [Fe/H] 
= -0.5 ± 0.09, [a/Fe] = 0.20 ± 0.10 for T eff = 4300 K. 
These results compare well with the adopted parameters of 
T efs = 5780 K, logg = 4.44, [Fe/H] = 0.0, [a/Fe] =0.0 for 
the Sun and T eff = 4286 K, log g = 1.66, [Fe/H] = -0.52 
and [a/Fe] = 0.33 for Arcturus (|Ramfrez fc Allende Prietol 
l201ll ). Note that the alpha enhancement values of Arc- 
turus in the literature show some scatter. When this pro- 
gram was done in 2010, we were aiming to match Arc - 
turus to the values of lPeterson. Dalle Ore, fc Kuruczl (|l993h 
of Teff = 4270 K, logg = 1 5, [F e/H] = -0.5 and [a/Fe] 
= 0.36. iNorris fc Da Costal (|l995h report an enhancement 
of 0.2 for the alph a elem ents of Ti, Si, Mg, closer to our 
derived value. iHilll (|l997l ) similarly report lower enhance- 
ment values. We could not reproduc e the higher value of 
IPeterson. Dalle Ore, fc Kuruczl ||l993h and similarly obtain 
a good match for our more alpha-enhanced open and globu- 
lar cluster calibrators. Our linelist is available in the online 
manuscript and a sample of the list is shown in Table [ 
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4.3 Radial Velocity Measurement 

Radial velocities were measured from the sky-subtracted and 
wavelength-calibrated spectra using the IRAF program fx- 
cor. For the template spectra, a restricted grid of the syn- 
thetic spectra was used. The fxcor internal errors for our 
spectra are typically 0.5 — 0.7 km s _1 , and the adopted tem- 
plate for each star was chosen to be the one for which the 
fxcor error was smallest. For fibre spectroscopy with a res- 
olution of 11,000 and mean S/N of 75, we would expect 
total radial velocity errors of about 1 km s _1 . To estimate 
our total errors empirically, we compared the weighted mean 
velocity and vel ocity dispersion of t he open cluster Melotte 
66 measured bv lSestito et al.l (|2008l ) (their radial velocities 
have errors of about 0.3 km s _1 ) with our measurement from 
10 stars in the cluster. The mean difference in velocity (Ses- 
tito - ARGOS) is -0.65 ± 0.39 km s" 1 , and the ARGOS 
velocity error for a single measurement is 0.92 ± 0.23 km 
s" 1 . 



4.4 Determining Stellar Parameters 

Our x 2 procedure to derive the stellar parameters follows 
Kirby et al. (2008) who measured the metallicity of red gi- 
ants from medium resolution spectra (R ~ 8000) over a wide 
metallicity range. This technique uses only weaker lines and 
was shown by Kirby et al. to reproduce the [Fe/H] mea- 
surements from high resolution spectra. Our spectra show 
weak lines of Fe I, Ti I, Si I, Al I, Mg I, O I and CN 
and the three strong Ca-triplet lines. Specific regions of 
the spectra, excluding the strong lines have been selected 
over which to perform the \ 2 minimization. Excluding the 
strong lines and using only the Fe and weak line regions 
means we return a true [Fe/H] measurement. Additionally, 
we cannot simultaneously match our strong lines over the 
wide parameter space of our calibration stars in our syn- 
thetic spectra so we exclude these fr om the analysis (see 
iKirbv. Guhathakurta. fe SnedenlbOOSf ). 

The procedure is iterative and requires at least four 
cycles to converge. The temperature derived from (J — K) a 
depends on the current estimate of metallicity and gravity; 
the gravity, [Fe/H] and [a/Fe] values derived from the \ 2 
process depend on the current temperature. The final cycle 
of the process gives the final adopted value of logg which 
comes from the correct a-enhanced grid. 

The a-enhancement level significantly affects the stel- 
lar line profiles for the a-elements, and it was essential to 
have grids of synthetic spectra that cover the full range of 
observed a-enhancement from to +1.0. 

Figure [6] compares observed and synthetic spectra for 
three stars from our fields including two near the extremes 
of our [Fe/H] distribution. The star in the lower panel has 
[Fe/H] = 0.3 and [a/Fe] = 0.15, the star in the intermediate 
panel has [Fe/H] = —0.4 and [a/Fe] = 0.3 and the star in the 
top panel has [Fe/H] ~ —2.5 and a high alpha enhancement 
[a/Fe] = 0.6. At this resolution and with a S/N = 50, the 
[Fe/H] value for the metal-poor star in the upper panel of 
Figure [6] is at the lower limit of measurability. For this star, 
with T e // ~ 5190 K, the Fe lines are barely visible against 
the noise, and the program cannot differentiate abundances 
lower than [Fe/H] = -2.5. For a S/N of 80, [Fe/H] values of 
~ —2.8 can still be measured. Stars of yet lower metallicity 



can be identified in our survey but their metallicities will be 
uncertain. 

We tested the dependence of the [a/Fe] ratio on S/N 
using synthetic spectra with noise added to the spectra and 
run through the \ 2 process. For stars with S/N = 50 per 
resolution element and [Fe/H] < —2.0, the uncertainty of the 
a-enhancement above our nominal reported error in [a/Fe] 
of 0.10 is increased by about 0.1 dex. For S/N = 80, it is 
increased by about 0.05. For stars with [Fe/H] > —1.0, there 
is no dependence of uncertainty on S/N for spectra with S/N 
> 50. 

The steps in the \ 2 process are: 

la) Choose wavelength regions for calculating \ 2 f° r deter- 
mination of [Fe/H] (weak Fe and weakly blended lines only). 
Exclude the Ca-triplet and lines of Ti, Al, Mg, Si. The re- 
gions used in the determination of [Fe/H] are shown in Fig- 
ure [6] There are about 30 Fe lines across this region. The 
equivalent widths of these lines is up to 500 mA for stars 
around solar metallicity and decreases to 200 mA for stars 
with [Fe/H] = -1.0. 

lb) Choose wavelength regions over which to calculate x 2 for 
measurement of [a/Fe] using a combination of the 17 avail- 
able weak Ti, Si, Mg lines in the synthetic spectra. These 
regions are shown in Figure [5] The equivalent width of the 
alpha-element lines is up to 500 mA at solar metallicity and 
about 200 mA at [Fe/H] = -1.0, similar to iron. 

2) Shift all observed spectra back to rest-frame wavelengths. 

3) Mask out any residual sky lines. 

4) Interpolate synthetic spectra on to the wavelength grid of 
the observed spectrum (0.24A channels) and measure reso- 
lution of observed spectra (R is nominally ~ 11,000 but 
varies by up to 10%). Then convolve synthetic spectra to 
resolution of the observed spectrum. 

5) Normalise observed and synthetic spe ctra in the same way 
using continuum regions identified from lHinkle et all (|2005l ) 
for Arcturus and the Sun. 

6) Make a first estimate of the temperature from the (J—K) 
colour. The interst ellar reddenin g for e ach star was de- 
termined from the ISchlegel et all (| 19981 1 reddening maps. 
The relat ionship between colo ur and temperature was taken 
from the lBessell et al] (| 19981 ) calibration, which was calcu- 
lated using the Castelli/Kurucz model atmospheres . Their 
calibration was transformed from the iBessell et al.l (|l998h 
(J — K) system to the 2MASS system usi ng the tranforma - 
tion (J — /sT)2mass = 0.983(J - .FObcp l| Carpenter! [2001 ). 
The Teff calibration depends on [Fe/H] and logg, and val- 
ues of [Fe/H] = —0.5 and logg = 2.5 were used to start the 
iteration. 

7) Compare each star with the synthetic spectra in the 
[a/Fe] = 0.0 grid and return the \ 2 value over the [Fe/H] 
regions from each comparison: 

2 \ ~* / Xi — Oi \ 2 

i=l 

where i — l...k are the wavelength channels for the [Fe/H] 
regions, Xi is the synthetic spectrum, Oi is the observed 
spectrum and Oi is the standard deviation of the observed 
spectra for a Poisson-noise dominated system. 

8) Interpolate to find the minimum value of x m t ne two 
dimensions of loggr and [Fe/H] to derive the best fit stellar 
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Figure 6. Ca-triplct region spectra for three stars in the sample. The star in the lower panel has T e ff= 4800 K, logg= 2.3, [Fe/H] 
= 0.3, [a/Fe] = 0.15. The star in the intermediate panel has T e f f = 4480 K, logg= 1.9, [Fe/H] = -0.4, [a/Fe] = 0.3 The star in the 
upper panel has T e ff= 5190 K, logg= 3.0, [Fe/H] = —2.5, [a/Fe] = 0.6. In each panel, the red spectrum is the synthetic spectrum for 
the adopted parameters and black dashed points is the spectrum as observed. The blue shaded regions correspond to the \ 2 regions 
to determine [Fe/H] and the purple shaded regions correspond to the regions to determine the [a/Fe]. Note the additional weak linos 
included near the Ca-triplet for stars < -1.75, shown in the top panel. The most metal-poor star in the top panel represents the limit of 
[Fc/H] measurement using the weak line regions. 



parameters. 

9) With t he new [Fe/H ] value for the star, redetermine T e ff 
using the iBessell et al.l |l998j) calibration. This second es- 
timate is the final temperature to be used, although the 
gravity will change slightly when the final log g value is de- 
termined from the \ 2 comparison of the star with the ap- 
propriate a-enhanced grid. 

10) Do a second x iteration over a narrower stellar pa- 
rameter range of the synthetic grid, around the minima of 
the first pass results, and interpolate over [Fe/H], log<? and 
[a/Fe], using the Fe regions for the [Fe/H] estimate and the 
a-regions of Ti, Si and Mg for [a/Fe]. 

11) . With this [a/Fe] value, use the appropriate a-enhanced 
grid to derive the final log g value. 

The procedure in step 11 was done in two independent 
ways which agree within ± 0.25 dex in log<?. The first im- 
plementation used the correct [a/Fe] grid and the whole 
spectral region to return the log g value. The reason for us- 
ing this method was to make use of the gravity information 
in the wings of the strong Ca-triplet lines. For this step, the 
log gf values for the triplet lines (included in the linelist) 
were adjusted so that the x 2 program returned the correct 
log<? values for the Sun and Arcturus for an [a/Fe] grid of 
and 0.2, respectively. The [M/H] returned from the x 2 com- 



parison to obtain the log g using the entire spectral region is 
disregarded as the \ogg and [M/H] could not be simultane- 
ously well determined over the wide parameter space of the 
calibration set. 

The second method fixed the [Fe/H] for the star and 
then interpolated only the x 2 iron and a-element regions 
shown in Figure [6] to determine the minimum x 2 for the 
single parameter logg. These weak line features are not as 
strongly gravity-sensitive as the Ca-triplet. However inter- 
polation in this single dimension enabled the logg value to 
be returned again to within ± 0.25 of the values of the first 
method. 

12) As described in section 4.1, the \ 2 program was checked 
against spectra of stars with known parameters over the 
abundance range [Fe/H] = —2.1 to 0.1. See the Appendix 
for comparison of the x 2 an d literature results for [Fe/H], 
[a/Fe] and logg for the set of calibration stars. 



5 ELIMINATION OF FOREGROUND DWARFS 

The low resolution spectra (R ~ 3000) from 5000 - 5600A 
were used for a first-pass elimination of foreground G and 
K dwarfs from the sample by computing an index that mea- 
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sures the strength of the gravity-sensitive Mgb/MgH fea- 
ture. The index was based on synthetic spectra of dwarf 
stars (log g ~ 4) and typical clump stars (log g ~ 2.7). How- 
ever, this was only partly successful. Due the large range in 
stellar parameters of stars in the sample, it was not possi- 
ble to cleanly discriminate against the dwarfs. In order to 
ensure that the more metal-rich ([Fe/H] > 0) giant stars 
were not excluded from the sample, it was necessary to set 
the index cut such that the more metal-poor dwarfs ([Fe/H] 
< —0.5) would also pass through. Additionally, for about 
30% of the stars, the signal level in the low resolution blue 
band was not high enough to allow secure discrimination 
of dwarfs and giants. The \ 2 program applied to the Ca- 
triplet region spectra was, however, easily able to identify 
the residual dwarfs in the sample. 



6 STELLAR DISTANCES 

Our selection criteria were aimed primarily at selecting red 
clump stars, because for these stars we can derive rela- 
tively accurate photometric distances. The sample does, 
however, include some fainter and brighter giants and also 
some foreground dwarfs. We use our estimates of stellar pa- 
rameters T e ff, logg, [Fe/ H] and [a/Fe] to identify the red 
clump stars. The study by Alvcs (2000) of a sample of Hip- 
parcos red clump stars with —0.6 < [Fe/H] < 0.1 gave a 
mean absolute K-magnitude of Mk = —1.61 ±0.03, and we 
adopt this mean value inde pendent of metallicity (see also 
iMcWilliam fc Zoccalil [2O10h ■ Alves fitted a gaussian to the 
luminosity function of his stars and found a dispersion in 
Mk of 0.22 mag. The mean error in the parallaxes is about 
5%, which contributes about 0.11 mag to the observed dis- 
persion of the red clump Mk values. For estimating the 
uncertainty in the distances of our clump stars, we adopt 
Alves's intrinsic dispersion of 0.22 mag for the Mk values . 
Thi s is similar to dispersions adop ted bv lStanek etafl (|l997f ) 
and lBabusiaux fc Gilmord (|2005l ). For the clump stars in our 
sample, the intrinsic dispersion of 0.22 in Mk is the domi- 
nant error source for the distance estimates. 

The mean logg for the red clump depends weakly on 
age and metallicity but is typically about 2.4. We can es- 
timate the expected intr insic spread i n log a from the dis- 
persion in Mk, using the lMarigo et al] (|2008T l isochrones for 
[Fe/H] = -0.2 and age 10 Gyr. Along the entire giant branch, 
the isochrones show that Mk increases almost linearly with 
logg, with dMjr/d(log g) = 2.3. The dispersion of 0.22 mag 
in Mk corresponds to a dispersion in logg of about 0.1. 
In the (logg, T e //) plane, the He-core burning red clump 
stars are superimposed on a background of first ascent gi- 
ants. This background of first giants can be minimised by 
choosing an appropriately narrow interval of logg: a fairly 
clean sample of red clump stars could be selected by tak- 
ing stars within say ±0.2 dex of the mean logg value, if the 
logg values were accurate enough to justify such a choice. 
Our measuring errors in logg are about 0.30 mag (see Ap- 
pendix), so we take stars with logg = 1.8 to 3.2 and T e // 
= 4500 to 5300 K as potential clump giants. 

In the above interval of logg and Te//, we have no way 
of telling which star is a red clump star and which is a first 
ascent giant. Near the mean logg of the red clump stars, the 
first ascent giants have almost the same value of Mk as the 
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Figure 7. The expected distribution of Mk for an unbiased pop- 
ulations of stars with [Fe/H] = —0.2 and age 10 Gy r lying on the 
giant branch with logg = 2 to 3 (derived from the lMarigo et aU 
(120081') isochrones). The luminosity function for the red clump 
stars has been convolved with a gaussian with a = 0.22 as for the 
Alves (2000) stud y. The mean Mk for the isochronc is slightly 
displaced from the lAlved feOOOh value. 



red clump stars themselves, but this is not so for first ascent 
giants with log g values towards the limits of o ur interval. As 
a guid e to the size of the effect, we can use the lMarigo et all 
(|2008l ) isochrones to estimate the probability distribution of 
Mk for the stars that lie on the giant branch and have 2 < 
log g < 3. This is shown in Figure[7]for a unbiased population 
with [Fe/H] = —0.2 and age 10 Gyr. Our sample is biased 
against the fainter stars because of the incompleteness of the 
2MASS survey at the fainter apparent magnitudes, so the 
numbers of intrinsically fainter giants with 2 < log g < 3 in 
our sample is reduced relative to the distribution shown in 
Figure [7] 

For the purposes of this survey, we take all of the poten- 
tial red clump stars to have Mk = —1.61 ±0.22, recognising 
that for some of the underlying first ascent giants the un- 
certainty will be somewhat larger. 

For stars that are not on the red clump, we derive Mk 
values using BaSTI isochrones of age 10 Gyr and appropri- 
ate [Fe/H]. The stars lie fairly close to the isochrones in the 
{T e ff, logg) plane (see paper III), and their absolute magni- 
tudes are derived by adjusting the temperature of the star 
to place it on the closest branch of the isochrone. The un- 
certainty in logg is 0.30 (see section 4.1) so the uncertainty 
in distance for non-clump giants is about 0.70 mag. 

In summary, for the clump star candidates which com- 
prise about 70% of our sample, the adopted distance error 
is ok = 0.22 mag (about 900 pc at the distance of the 
Galactic centre). For the non-clump giants, the distance er- 
ror is ok = 0.70 mag, dominated by the uncertainty in the 
logg estimates. The near-linear relationship between Mk 
and log g on the giant branch means that the error distribu- 
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tions for the distance moduli for both clump and non-clump 
giants are close to Gaussian. 



7 CONCLUSION 

In this paper we have presented the ARGOS survey. Pa- 
pers III (metallicity components) and IV (kinematics) will 
show that we have achieved our goals of finding the red 
clump giants in the bulge and discovering the small frac- 
tion of metal-poor stars in the bulge region: ~ 0.15% of our 
stars have [Fe/H] < —2.0. These later papers will describe 
the components of the metallicity distribution function that 
are seen throughout the bulge region, discuss their origin in 
the context of bulge formation theory, and show how the 
kinematics of the bulge stars change with metallicity. 
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APPENDIX A: APPENDIX 

The uncertainties in our estimates of [Fe/H], [a/Fe] and logg 
were determined by comparing our parameter estimates with 
those in the literature for the Sun, Arcturus, 16 bulge stars 
observed by Zoccali et al (2008), the open cluster Melotte 
66 and six globular clusters. Figures A1-A3 show the dif- 
ferences A between our \ 2 estimates and the literature esti- 
mates for [Fe/H] and [a/Fe], and between our \ 2 an d com- 
parison estimates for logg. For the [a/Fe] ratio we use the 
mean of the Mg, Si and Ti abundances. 

The scatter of the differences A comes from our \ 2 mea- 
suring errors and errors in the comparison values of the pa- 
rameters 

2 2 2 

C"A = 0" x 2 + O" comp 

where cta is the weighted dispersion of the differences A. 
For [Fe/H], cta = 0.13 and we estimate that a comv — 0.10, 
so a x 2 = 0.09. For [a/Fe], we find that a x 2 = 0.10. 

For the log g errors, some of our comparison stars are in 
clusters and we compared o ur y 2 estimates with log g values 
from isochrone fits using the lMarigo et al] (|2008l ) isochrones. 
Colour and magnitude data and the adopted cluster ages 
and distances were taken from the literature. For example, 
the age of Melotte 66 was taken to be 4 ± 1 Gy r and its dis- 
tance modulus (m — M) = 13.2 (+0.3, —0.1) |Kassis et al.1 
1 19971 ). The stars fit well on the isochrone after small adjust- 
ments within the errors were made to the distance. The main 
source of error in the log g determination for this procedure 
is the uncertainty in the age of the population. For Melotte 
66, an age uncertainty of ±1 Gyr corresponds to an uncer- 
tainty <7io g9 = 0.20 dex. For the older clusters, an age error 
of 1 Gyr corresponds to a lower uncertainty <7i og9 = 0.13. 
From the distribution of differences shown in Figure A3, we 
estimate that our uncertainty on logg from the \ 2 process 
is criogg = 0.30 dex. 



3 There are two lZoccali et all j200Sf) stars that lie clearly above 
the 1 — <j crrorbars of the distribution in the figure. One of these 
stars falls on the RGB br anch from the ARGO S calibration and 
the AGB branch for the IZoccali et al. calibration. The 

oth er outlier is a poor m atch to the isochrone for both ARGOS 
and lZoccali et at calibrations. 
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Melotte 66, [Fe/H] = -0.3, S/N = 60 
47 Tuc, [Fe/H] = -0.75 
NGC288, [Fe/H] = -1.3, S/N = 33 
NGC1904,[Fe/H] = -1.6, S/N = 34 
NGC 2299, [Fe/H] = -1.9, S/N = 55 

Arcturus 



i 



<B -0.1 

"3 



0.2 0.3 0.4 

[a/Fe] from literature 



Figure Al. [a/Fe] calibration. For the open and globular cluster data, the error bars in the x-direction represent the measurement 
error of the literature comparison values <T com p, and the error bars in the y-direction represent the combined internal x 2 error and the 
measuring errors for the comparison values, cr? = {<J 2 omp + c 2 2 ). For Arcturus and the Sun, the errors in the y-direction are the internal 
X 2 fitting errors. 




Figure A2. [Fe/H] calibration. The error bars in the x-direction represent the estimated literature errors. As for Figure Al, the error 
bars in the direction of the y-axis represent the combined literature error and the x 2 fitting errors for the individual stars and the 
measuring error for the open and globular clusters, a 2, of their mean [Fe/H]. 
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* Melotte66, ilx'/HJ ■ 0.3, S/N 60 
» 47 Tuc, [Fe/H] = -0.75, S/N = 80 

• Zoccali stars, [Fe/H] - -1.2 to 0.2, S/N = 77 

■ NGC288, [Fe/H] - -1.3, S/N - 33 
A NGC1904,[Fe/H] = -1.6, S/N = 34 

■ NGC 2298, [Fe/H] = -1,9, S/N = 55 
□ Sun 

♦ Arctums 



3.5 4.0 4.5 

log g from isochrone/literature 

Figure A3, logg calibration. The error bars in the x-direction represent the uncertainty in the estimated logg for the calibration stars. 
The error bars in the y-dircction represent the combined error of the x 2 program and the uncertainty of the log g value for the comparison 
obj ects. Note the atypic ally large errors from the poor S/N spectra for the cluster NGC1904. The logg values for the Sun and Arcturus 
and lZoccali et alj {2008h are taken from the literature and all logg values for cluster stars are from isochrones estimates (see text) 3 . 

Table Al. Comparison of stellar parameters for the 16 stars observed for ARGOS measured bv lZoccali et all | j2008t ). 







Zoccali et al. (2008) 






ARGOS 




OGLE 


S/N 


T eff 




[Fe/H] 


T eff 


log 9 


[Fe/H] 


[cc/Fe] 


100047c6 


80 


4350 


1.72 


-0.33 


4470 


1.7 


-0.26 


0.3 


104943c6 


75 


5300 


2.1 


-1.2 


5660 


3.4 


-0.90 


0.2 


75097c7 


65 


4500 


1.8 


-0.09 


4650 


2.2 


0.10 


0.1 


57883c7 


65 


4550 


1.8 


-0.25 


4440 


1.6 


-0.44 


0.3 


90337c7 


70 


4850 


2.02 


-0.27 


4770 


2.3 


-0.36 


0.2 


64860c7 


75 


4500 


1.8 


0.09 


4790 


2.2 


0.26 


0.2 


98458c7 


80 


4300 


1.68 


-0.02 


4540 


2.4 


0.1 


0.1 


87242c8 


70 


4300 


1.67 





4440 


2.0 


0.1 


0.2 


43791c3 


65 


4950 


2.2 





5160 


2.6 


0.0 


0.1 


32080c3 


85 


5000 


2.0 


-0.17 


5450 


3.1 


-0.12 


0.2 


23017c3 


80 


4250 


1.65 


0.09 


4410 


1.6 


0.18 


0.1 


38354c3 


75 


4700 


1.83 


-0.61 


4620 


1.9 


-0.68 


0.2 


76187c3 


80 


4550 


1.79 


-0.42 


4700 


2.1 


-0.34 


0.4 


56410c2 


75 


4600 


1.81 


-1.1 


4570 


1.3 


-1.20 


0.5 


19402cl 


60 


4550 


1.83 


-0.61 


4730 


2.3 


-0.40 


0.2 


212324c6 


75 


4800 


1.88 


-0.32 


4960 


2.4 


-0.28 


0.2 



Table A2. Sample of the linelist available in full in the electronic version of this paper. Following iKirbv. Guhath akurta. &; Sn cdcn 
j200Sh . changes from original log gf values are indicated by values with single decimal places. The linelist is in MOOG readable format. 



Wavelength (A) 


Species 


EP (eV) 


log gf 


8467.937 


20.000 


4.554 


-3.313 


8468.407 


26.000 


2.223 


-2.3 


8468.470 


22.000 


1.887 


-1.2 


8468.592 


26.100 


9.700 


-3.0 


8468.657 


607.00 


0.916 


-1.597 


8468.732 


28.100 


6.952 


-4.8 


8468.832 


25.000 


6.412 


-2.8 


8468.849 


12.000 


5.932 


-2.240 


8468.872 


607.00 


0.916 


-2.9 
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